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ABSTRACT 

We present a detailed analysis, in the framework of the MFA approach, of 
the critical behaviour of the lattice Schwinger model with Wilson fermions on 
lattices up to 24^, through the study of the Lee- Yang zeros and the specific 
heat. We find compelling evidence for a critical line ending at k = 0.25 
at large (3. Finite size scaling analysis on lattices 8^,12^,16^,20^ and 24^ 
indicates a continuous transition. The hyperscaling relation is verified in the 
explored /3 region. 
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1. Introduction 

Although apparently far from physical reality, the Schwinger model, i.e., 
the theory describing the interaction between photons and electrons in 1 + 1 
dimensions, has been a favourite framework for theoretical and numerical 
exercises for at least two reasons. 

First, the massless Schwinger model can be analytically solved in the 
continuum and this is the reason why it has always been used as a laboratory 
for the development of numerical algorithms for dynamical fermions. 

Second, and due to the special low dimensionality dynamics, this model 
shares many physical properties with four dimensional QCD, the gauge the- 
ory describing the strong interaction of hadrons. In fact the charge is confined 
in the Schwinger model and long range forces are absent in it. Of course the 
Schwinger model has also some characteristic features which are not common 
to QCD like the property of superrenormalizability or ultraconfinement (it 
effectively describes free bosons) and the absence of the Goldstone boson. 
However we expect a phase diagram in the lattice regularized Schwinger 
model qualitatively similar to that of lattice QCD, in particular since, as 
can be derived from dimensional arguments, the continuum limit is reached 
at infinite gauge coupling. This is one of the motivations for this paper. 

The Schwinger model on the lattice, using the Kogut-Susskind regular- 
ization for the fermion fields, has been the subject of extensive analysis. The 
phase structure of the model in the gauge coupling (3 and fermion mass m 
plane is well known and the continuum value for the chiral condensate has 
been reproduced within three decimal places , the best numerical result to 
our knowledge, using the MFA approach 0. In the Wilson regularization 
for fermion fields however, the situation is not so clear. In this scheme, chiral 
symmetry is explicitly broken even for vanishing bare fermion mass and this 
is the unavoidable price to pay to overcome the species doubling problem. 
Therefore no order parameter can be used for analyzing the phase diagram, 
which besides the larger number of degrees of freedom if compared with the 
Kogut-Susskind regularization, makes the analysis far from trivial. 

Furthermore almost all the standard algorithms to simulate numerically 
dynamical fermions, like hybrid Monte Carlo, which are based on the inclu- 
sion of the fermion determinant in the integration measure, need to dupli- 
cate the number of fermion species in order to avoid negative values for the 
determinant of the Dirac operator, which makes impossible to investigate 
the phase structure of the one flavour Schwinger model. This is the reason 
why until very recent time, when an analysis of the phase diagram of the 
Schwinger model in the Wilson regularization using non standard methods 
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was performed 0, very little was known about it. 

In this article we will report the results of a numerical simulation of the 
one flavour Schwinger model in the Wilson regularization, by means of MFA 
simulations. The structure of the Lee- Yang zeroes in the complex hopping 
parameter k plane strongly suggests the existence of a critical line: in our 
simulations the critical k is 0.27 at = 2, ending at k = 0.25, (3 — oo, the 
critical point where the continuum limit is recovered. We also describe results 
concerning the specific heat showing a sharp peak along a line Kc = fiPc) 
in the k, (3 plane, in good numerical agreement with predictions from the 
Lee- Yang zeros analysis. The results of the finite size scaling analysis for 
the Lee- Yang zeros and chiral susceptibility allow the determination of the 
critical exponents, and the fulfillment of the hyperscaling relation gives a 
definite confirmation that we are dealing with a real second order phase 
transition line. 

2. Theoretical Grounds 

The lattice action for the massive Schwinger model with Wilson fermions 
is given by 

S{(3,k) = Sf{k)+(3Sg (1) 

where Sf{k) contains the kinetic and mass terms for the fermion field as well 
as the fermion-gauge interaction term, and we have chosen for the pure gauge 
action Sq the standard noncompact regularization for the abelian model. The 
fermionic action Sf{k,) reads 

fl,X 

+tP{x){l--f^)U^{x)^{x + fi)] - J24j{x)^{x) (2) 

X 

Equation (2) defines the Dirac matrix operator A which can be written 

as 

A = -I + kM (3) 

where / is the unit matrix and M a matrix of dimension equal to the lattice 
volume times the number of Dirac components. 

The use of the MFA approach in the numerical simulations of the system 
described by action (1) has at least two very important advantages when 
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compared with other standard approaches. First, since MFA is based on 
the computation of the fermion effective action defined as the logarithm of 
the mean value of the fermion determinant at fixed pure gauge energy and 
this mean value is positive definite at least in the physically interesting region, 
we have no problems to simulate the one flavour model. Secondly, since the 
K dependence factorizes in the non trivial part of the fermionic operator (3), 
we do not need to repeat the numerical simulations when changing the value 
of the hopping parameter k, which is mandatory for exploring the zeros of 
the partition function in the complex k plane. Let us notice also that, while 
MFA has been extensively used and checked in lattice gauge theories with 
Kogut-Susskind fermions, this is the first time it is applied to a lattice model 
with Wilson fermions. 

The technical details of MFA can be found by the interested reader in 
refs. 0. As stated before, all the applications of MFA to lattice gauge 
theories in these references were done using staggered fermions. The main 
difference between staggered and Wilson fermions comes from the different 
structure that the Dirac operator has in each formulation. In the Wilson 
case, where the fermion matrix has the structure given by equation (3), the 
MFA approach works as follows: first we generate well decorrelated gauge 
field configurations at fixed noncompact gauge energy 

E = ^ (M^) + + f^)- + ^)- M^)) (4) 

and then we diagonalize the M matrix for each generated configuration. Due 
to the non hermiticity of the matrix we can not use the standard Lanczos 
algorithm; the eigenvalues are found using a standard library diagonalization 
routine. From the eigenvalues of the M matrix we can reconstruct the de- 
terminant of the Dirac operator A at any value of the hopping parameter k 
in a trivial way. The partition function associated to the action (1) can be 
written then as 

Z{P, «:) = ^ = / t/En(E)e-^^^ ^ (5) 

where n{E) is the density of states at fixed pure gauge energy and Cn{E) 
in (5) stands for the mean value of the nth coefficient of the polynomial 
describing the fermionic determinant, the average being computed over gauge 
field configurations at fixed pure gauge energy, i.e.. 
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Cn{E) = [dA^]Cn{A^{x))5{ {A^{x)+A,{x+fi)-A^{x+u)-A,{x))-VE) 

(6) 

Since the density of states at fixed pure gauge energy 

n{E) = I [dA^]6{ J2 (M^) + + /^) - + M^)) - VE) (7) 

can be analytically computed in the noncompact model, we can reconstruct 
the partition function (5) from the knowledge of the coefficients Cn{E) through 
interpolation and one-dimensional integration. Had we used the more stan- 
dard compact Wilson formulation for the gauge fields, the standard procedure 
would be very similar with the only difference that in such a case, the density 
of states, in general, would have to be computed numerically. 

If on the other hand we are interested in the computation of vacuum ex- 
pectation values of physical operators like the chiral condensate, the standard 
procedure is the one described in 0. We will not repeat here the details of 
this procedure but only will remember that it is based on the computation of 
the fermion effective action as a function of the gauge energy and the com- 
putation of the mean value of the operator times the fermion determinant 
over gauge field configurations of fixed pure gauge energy. 

3. The Lee- Yang zeros in the complex k plane. 

The first step to the determination of the position of the Lee- Yang zeros 
in the complex k plane is the computation of the coefficients Cn{E) of the 
averaged determinant of the Dirac operator. For the calculation of Cn{E) 
we proceed as follows: first we chose a set of values of energy, in the range 
selected to cover the support of the weight function in (5) for the values of 
f3 we are interested in. Then for every value of E in the set we generate 
gauge field configurations using a microcanonical code; the generation of 
gauge fields at fixed energy is not the costly part of the whole procedure, so 
we can well decorrelate the configurations used for measuring the fermionic 
operator. Then, as stated before, we compute exactly the eigenvalues of the 
M matrix from which we reconstruct by a recursion formula the coefficients 
of the fermionic determinant. 

At the end we have the coefficients Cn{E) evaluated at discrete energy 
values: a polynomial interpolation allows the reconstruction at arbitrary 
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values of the energy E, in order to perform the numerical integration in (5) 
and obtain the coefficients C„(/?) that can be regarded as the final product 
of this part of the numerical procedure for the determination of the Lee Yang 
zeros. These coefficients are then used for the determination of the roots of 
the polynomial Z{f3, k). 

The main features of Z{(3, k) are the following: in a volume V, the par- 
tition function is a polynomial of order N = IV in k and the typical range 
of the coefficients is of order e^. Thus, the main numerical problem is the 
efficiency of standard root finders in the determination of the zeros. Here we 
have used a method developed in order to analyse the partition function 
zeros in the four dimensional compact model. This algorithm is based 
on well known properties of analytic functions on the complex plane which, 
in particular, allows the determination of the number of zeros for a given 
function inside a region of the complex plane, provided the function has no 
singularities inside this region. 

In Figure 1 we plot the location of all the zeros (in the complex k plane) 
in a 16^ lattice at /3 = 7. Figures 2 a,b contain the zeros closest to the real 
axis in all the lattices used at /3 = 7, 10 

To estimate the statistical errors on the position of the Lee- Yang zeros we 
followed the procedure described in [Q: a standard jack-knife method is used 
to produce n averaged partition functions and n estimates of the location 
of a given zero. The largest error on the distance of critical zeros (those 
with smallest imaginary part) with respect to the free fermion critical point 
(0.25,0.0) is of order 2%. In all the cases reported in these figures, we can 
see how the nearest zero to the real axis approaches it with increasing lattice 
size, thus suggesting the existence of a phase transition line in the (3 — k 
plane. 

In Figure 3 we show the phase diagram in the (3 — k plane obtained from 
8^, 16^ and 24^ lattices. The real part of the location of the zero lying nearest 
to the real axis gives an estimate of the value of the critical k. Note that the 
critical line moves upward with increasing volume. 

A detailed analysis of the scaling behavior for the Lee- Yang zeros and Kc 
with the lattice size will be presented in section 5. 

At smaller values of (3 our procedure has some problems: the signature 
is the appearence of unphysical real zeros at finite volume, so that the es- 
timation of the critical point becomes unsafe. This behaviour is related to 
the large fiuctuations in the averaging procedure for the coefficients, already 
pointed out in 0. Nevertheless, if we take the zeros with smaller, but not 
zero, imaginary part as the critical zeros we can continue the critical line 



6 



down to P = 0, ending at kc = 0.36 for a 8^ lattice and kc = 0.35 for a 16^ 
lattice. Note that in 0] an exact result for the partition function at /3 = 
for a 8^ lattice is reported, giving a critical point at kc = 0.377. 



4. The specific heat. 

Chiral symmetry is always explicitly broken and there is no order parame- 
ter in this realization of the model, so n cannot be identified with an external 
field. We take it instead as a temperature, defining the the associated specific 
heat as usual: 

where T is the free energy density. 

The specific heat is related to the chiral susceptibility defined as 

X = -2k^-^ < ijj(x)-ip(x) > (8) 
an 

differing only through a regular function of k which does not infiuences the 
critical behavior. The susceptibility, which is what we have measured, di- 
verges at the transition in the thermodynamical limit, the divergence becom- 
ing a sharp peak on finite lattices. 

The study of the position of the maximum at different lattice sizes gives a 
way to search for the phase transition line which is numerically independent 
from that presented in the previous section. We remind that, as the model 
is not chirally invariant, the analogy with a magnetic system fails. 

Figures 4a,b contain the results for the susceptibility (8) at two typical 
values of the gauge coupling j3 against k in 8^, 12^, 16^, 20^ and 24^ lattices. 
As expected for a real phase transition, the susceptibility shows a well defined 
maximum at some critical value k^, the height of these peaks increasing with 
the lattice size and eventually diverging in the infinite volume limit. If we 
define the critical Kc at each f3 as the value at which the susceptibility x 
takes its maximum value, we get for the lattices analyzed the phase diagram 
reported in Figure 5, in very good agreement, in the largest lattices, with the 
one obtained from the analysis of the Lee- Yang zeros. 



5. Finite size scaling analysis. 

We will now present a detailed analysis of the scaling behavior for the 
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Lee- Yang zeros and specific heat, mainly based on simulations on 16^, 20^ 
and 24^ lattices. 

In the analysis of the phase diagram from the location of the Lee- Yang 
zeros of Figure 3, the real part of the zero lying next to the real axis defines 
the critical k at each lattice size. In order to have a real phase transition, the 
imaginary part of the critical zero should vanish in the infinite volume limit. 
This has been explicitly checked by assuming that the imaginary part of the 
critical zero Zc{L) as a function of lattice size is described by the function 

Tm Zc{L) = ao + aiL~~ . (9) 

In all the f5 region explored, the value of Oq is compatible with zero. The 
scaling analysis at this point is made using the relations 

Im Zc{L) L-^ . (10) 

for the real and imaginary parts of the critical zero. 

In Figure 6 we show the typical scaling behavior as a function of the 
volume for the imaginary part of the critical zero, at /5 = 10. Similar results 
have been obtained at different (3. From the slopes of the lines fitting the 
data from 16^, 20^, 24^ we determine the scaling exponent z/ as a function of 
(3. This is shown in Figure 7. The value of the critical index clearly indicates 
that the transition is continuous at any value of f3. We note, moreover, that 
our results are compatible with z/ = |, at any /3. 

Concerning the susceptibility the critical k reported in Figure 5 was 
defined, as stated before, as the value of k at which the susceptibility takes 
its maximum value. Both the height H{L) of the peak as well as its position 
Kc{L) depend on the lattice size L. Standard finite size scaling theory tell us 
that these quantities scale with the lattice size as 

H{L)r^0 (11) 

where a in (11) is the specific heat exponent. 

In Figure 7 we present also the scaling index of as a function of f3. 
Its value is again consistent with a continuous phase transition. The value 
obtained implies a ~ z/. 

Since in the critical region the singular part of the free energy behaves 

as: 
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Ts^n9 ~ («: - «:c)'-" (12) 

and given that the only relevant length in this region should be the correlation 
length, we get 

T,,ns ~ r" ~ - /^c)'" (13) 

from which the standard hyperscaling relation between the correlation length 
and specific heat exponents follows 

dv = 1-a (14) 

Our data for the critical exponents satisfy relation (14) within the errors, 
in the whole 5 < /? < 12.5 region; the typical deviation from (14) is of the 
order of 0.07. 

We can conclude from the finite size scaling analysis of the Lee- Yang 
zeros and chiral susceptibility that the lattice Schwinger model with Wilson 
fermions has a real continuous phase transition ending at /3c = C)0, Kc = 0.25. 
We are able to derive the critical exponents of the model and we find perfect 
agreement with the hyperscaling relation. 

In Figure 8 we show the scaling behavior of the real part of the critical 
zero at (3 = 10. The infinite volume value of Kc can be inferred using the 
Finite Size Scaling in (10). In principle this analysis, carried out for several 
values of /3, could give an estimate of the infinite volume critical 
function of the gauge coupling constant. However it seems very difficult to 
get reliable values for Kc(oo) from the relation (10) for the real part of the 
critical zeros since this is not a universal relation. 

6. Discussion. 

In the previous sections we have shown that the Schwinger model regular- 
ized on a lattice and with Wilson fermions has a continuous phase transition 
in the {(3, k) plane, where the correlation length diverges. The finite size 
scaling analysis shows a correlation length exponent u taking a value around 
2/3 along this transition line, a result which appears very reliable since, as 
previously shown, the hyperscaling relation (14) is always verified. This value 
of u is in contrast with the value u = 1 obtained at the end point /5 = oo of 
the transition line as well as with the value, again z/ = 1, at /5 = suggested 
by the analysis reported in 10] We would like to notice that violations 
to Universality in fermionic systems, probably due to the long range forces 
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induced by the fermion fields, have been previously observed in the Gauged 
Nambu-Jona Lasinio model 0, [10| as well as in a fermion-gauge-scalar model 



Let us finally comment on the physical meaning of the continuous phase 
transition line. In QCD4 this transition line, whose existence and location 
are not as clear as here, is assumed to be the line along which the pion is 
massless. In two-dimensional models, as well known, Goldstone bosons are 
absent. Moreover, the correlation length which diverges in our case is that 
associated to the scalar particle, not the pseudoscalar one. Nevertheless, it is 
our prejudice that this is the line along which the fermion remains massless 
even if, due to the pathologies of two-dimensional models, massless scalars 
or pseudo-scalars are absent. 
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Figure captions 



Figure 1. Position of all the zeros at /5 = 7.0, 16^ lattice. 

Figure 2. Zeros closest to real axis, /3 = 7, 8^, 12^, 16^, 20^, 24^/3 = 10 

(b). 

Figure Z. (3 — k phase diagram from Lee- Yang zeros. 

Figure 4. Chiral susceptibility versus k (3 — 6.9 (a) and /3 = 9.9 (b); the 
error shown is the tipical one. 

Figure 5. {3 — k phase diagram from chiral susceptibility (upper curves), 
8^, 16^, 24^ lattices, the lines are the data of Figure 3. 
Figure 6. Imaginary part of the critical zero versus lattice size L, (3 = 10.0. 
Figure 7. i/ , ^ versus /3, 16^, 20^, 24^ lattices. 

Figure 8. Real part of the critical zero versus L,P — 10.0; the hue is the 
best fit using (10). 
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